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ABSTRACT 

The fundamental plane of early-type galaxies is a rather tight three-parameter correlation 
^J ' discovered more than twenty years ago. It has resisted a both global and precise physical 

C) \ interpretation despite a consequent number of works, observational, theoretical or using nu- 

merical simulations. It appears that its precise properties depend on the population of galaxies 
in study. Instead of selecting a priori these populations, we propose to objectively construct 
homologous populations from multivariate analyses. We have undertaken multivariate cluster 
and cladistic analyses of a sample of 56 low-redshift galaxy clusters containing 699 early-type 
galaxies, using four parameters: effective radius, velocity dispersion, surface brightness aver- 
aged over effective radius, and Mg2 index. All our analyses are consistent with seven groups 
that define separate regions on the global fundamental plane, not across its thickness. In fact, 
each group shows its own fundamental plane, which is more loosely defined for less diversi- 
fied groups. We conclude that the global fundamental plane is not a bent surface, but made of 
I> ' a collection of several groups characterizing several fundamental planes with different thick- 

••/"") , nesses and orientations in the parameter space. Our diversification scenario probably indicates 

that the level of diversity is linked to the number and the nature of transforming events and 
that the fundamental plane is the result of several transforming events. We also show that 
our classification, not the fundamental planes, is universal within our redshift range (0.007 - 
lO ' 0.053). We find that the three groups with the thinnest fundamental planes presumably formed 

through dissipative (wet) mergers. In one of them, this(ese) merger(s) must have been quite 
ancient because of the relatively low metallicity of its galaxies, Two of these groups have 
subsequently undergone dry mergers to increase their masses. In the k-space, the third one 
clearly occupies the region where bulges (of lenticular or spiral galaxies) lie and might also 
have formed through minor mergers and accretions. The two least diversified groups probably 
did not form by major mergers and must have been strongly affected by interactions, some of 
the gas in the objects of one of these groups having possibly been swept out. The interpre- 
tation, based on specific assembly histories of galaxies of our seven groups, shows that they 
are truly homologous. They were obtained directly from several observables, thus indepen- 
dently of any a priori classification. The diversification scenario relating these groups does 
not depend on models or numerical simulations, but is objectively provided by the cladistic 
analysis. Consequently, our classification is more easily compared to models and numerical 
simulations, and our work can be readily repeated with additional observables. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: fundamental parameters - methods: statistical 



1 INTRODUCTION and on the other hand between their effective rad ius R e and surfac e 

brightness averaged over effective radius < fi e > JKormendvll 19771) . 
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Physical understanding of astrophysical objects most often uses 

correlation diagrams. For early-type galaxies, such scaling laws . . . . , . , .. , „ , 

- ■ ' - ing a three-parameter relation of the form : log it,, = a logo - + b < 



These correlations are rather tight, but the scatter is still reduced us- 



. J J1 ° , ina a three-parameter relation ol the torm : log K e = a loser + b < 

have been for instance established on one hand between optical lu- " hr : — n|,„„j itt , . „ _ — nnTCS^k m - 

. . , t , , . ,. J c u „ t i l im^ V' > + c toressleretal]|l98llDiorgovski&Davislll987h. This 

minosity and central velocity dispersion <r(Faber & Jackson 1976), . , „ . ^ , , , . . ... ' n 

relation extends to taint and low-mass galaxies (e.g. INieto et all 
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1990). This is the fundamental plane (hereafter FP) of early-type 
galaxies. 
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A long-running difficulty is the so-called "tilt" of the FP with 
respect to the "virial plane" obtained with the virial theorem and 
some simple assumptions about the population of early-type galax- 
ies. Indeed, thi s tilt is different for di f ferent types of galaxies, like 
disk ones (e.g. iRobertson et alj|2006l : iHopkins et alj [2008). Many 
studies have been devoted to this problem without a clear solution. 
The motivation is to obtain pure correlations both to constrain the 
models better and to use them as probes of characteristics which are 
difficult to measure or strongly biased. In particular, the FP could 
in principle be a powerful tool to measure distances. But a proper 
calibration is required, and this appears difficult with its different 
tilts depending on the galaxy populations. 

Interpretation of th e FP very often assu mes some homol- 
ogy which is defined by iGargiulo et al.l ( 120091) as: "systems with 
density, luminosity and kinematics structures equal over the en- 
tire early-type sequence and with constant mass-to-light ratios". 
van Do kkum & Stanford (2003) define homology in such a way 
that the evolution of the FP is due only to the evolution of M/L. 
These definitions are certainly linked to the assumption that all 
early-type gala xies are assembled in the same way, like dissipa- 
tional mergers JRobertson et all 2006; Hopkins et al. 2008). 

Reality is however more complicated (see e.g . 
Bender et alj Il992l; Ijorgensen et alj 1 19961 : iBorriello et all 120031 : 
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very probably not all the result of the same formation process. 
Basically, what is needed is some invariant that allows us to 
trace a given object or class of objects through changes due to 
evolution. This invariant has been hoped to be the FP relation with 
some universality among a given population, universality which 
is implied by the definition of homology above and characterized 
essentially by M/L. This provides a rough way of simplifying the 
many variables that can evolve and may hide important sides of 
galaxy evolution. 

Models assuming only one homologous population have ap- 
parently failed to fit the FP in its entire extent. Selection criteria 
were proposed in order to obtain more homogeneous samples and 
thus define what could be called a "purer" fundamental plane. How- 
ever, these criteria are necessarily arbitrary, subjective and/or model 
dependent. The difficulty is that many parameters are known to in- 
fluenc e the global shape of the FP. Tilt of the FP ( IRobertson et alj 
2006), warps, dispersion, c hanges with redshift, depen dence of the 
mass-distribution on mass JNigoche-Netro et aU 2009) among oth- 
ers, show that even though the FP looks tight in the log cr, < p. e > 
and \o%R e space, additional parameters could still play a role. 

The universality of the FP is also questioned: is it a bent plane 
( IGargiulo et al.l 20090 or a bent surfac e approximated by a collec- 
tion of planes JP'Onofrio et al. 2008)? Obviously, this problem is 
related to the choice of the sample, that is to the definition of ho- 
mology. Theorists cope with many parameters that may influence 
the evol ution of a given galaxy , especially when mergers are con- 
sidered JRobertson et al .1120061) . Testing parameters one after the 
other, both theoretically and observationally, takes a lot of time, 
and might partly explain why after so many years, this tight corre- 
lation keeps most of its mysteries. 

We think that it is time to explore new methodologies to better 
characterize and understand the FP relation. Multivariate clustering 
approaches are more objective in selecting really multivariate "ho- 
mologous" sub-populations of galaxies. This requires to explicitly 
assume that the "fundamental plane" is a priori not universal, and 
to understand it as a correlation in the (log R e , log cr, < fj e >) space 
that could depend on the sub-population. Since galaxies are evolu- 
tive objects, homology can be more rigorously defined by 'similar- 



ity due to same class of progenitor' . It is the astrophysical equiv- 
alent of 'similarity by common ancestry' in cladistics, a statisti- 
cal method designed to relate evolutionary objects and developed 
mainly by biologists. The use of many parameters is necessary to 
find the true homology, and to prevent analogy (same characteris- 
tics due to convergent evolution) to yield false lineages of galaxies. 
Cladistics does not assume a priori properties linked to homology, 
but rather relies on all the pertinent parameters to construct homol- 
ogous groups. 

In the present paper, we have performed multivariate clas- 
sifications with two independent approaches, cluster analysis and 
cladistic analysis, of a carefully chosen sample from the literature. 
For this first kind of study, we consider only the three parameters 
of the FP, log cr, < /i e >, log R e , plus Mg2 that are all given for this 
homogeneous sample. The first justification is that the number of 
parameters is small, so it is not necessary to use PCA because we 
can physically discuss the variation in the light of these few param- 
eters easily. The second advantage is that we are interested to study 
the evolution of early-type galaxies in the light of their FP proper- 
ties, and most of the authors have discussed evolution of galaxies 
on the basis of global FP. 

The first approach we have used is by a multivariate 
technique known as K-means Cluster Analysis using the pa- 
rameters above. The second approach, known as astrocladis- 
tics, is based o n the evolutionary nature of both galaxies and 
their properties l[F raix-Burnet et al. 2006a bug: IFraix-Burnetll2009l : 
iFraix-Burnet et alj|2009r) . The clustering technique compares ob- 
jects for their global similarities, while astrocladistics gathers ob- 
jects according to their "histories". The two techniques are indeed 
complementary, the first one identifying coherent groups, the sec- 
ond one establishing an evolutionary scenario among groups of ob- 
jects. They are also totally independent, so that the comparison of 
their results is extremely instructive from a statistical point of view. 

This paper is organized as follows. The data and the different 
multivariate analyses we performed are described in Sect. [2] Com- 
parison of the structures resulting from these analyses, character- 
ization of the fundamental planes of individual groups, evolution 
properties within the global FP and properties of the groups, are 
presented in Sect. [3] The discussion on the meaning of the groups 
we have found is given in Sect. [4] before our conclusion in Sect. [5] 
We present detailed descriptions of the cluster and cladistic analy- 
ses respectively in Appendix lAl and Appendix|B] and provide addi- 
tional diagrams in Appendix Icl 



2 MULTIVARIATE ANALYSES 

2.1 Data 

For this first study, exploratory, it was important to choose the sam- 
ple carefully to be able to understand the result. We first need a 
homogeneous sample, devoid as much as possible of systematic bi- 
ases. We need to remain at low redshifts, but ideally integrating 
several clusters. We are limited in the number of objects in the as- 
trocladistic analysis. At the current stage of development of this 
novel approach, the reasonable limit is a thousand objects, and it is 
clear that a sample with 10 5 objects is untractable directly. 

We thus chos e the data compiled and standardized by 
iHudson et alj J200ll) for 699 early-type galaxies in 56 clusters, in 
the redshift range from 0.007 to 0.053 (SMAC catalog). This sam- 
ple has been carefully established in order to obtain distance esti- 
mates from the FP and study streaming motions in galaxy clusters 
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( Hudson et alj|2004l) . The data include the three parameters of the 
FP space: R e in arcsec, < p e > in the R-band in mag arcsec~ 2 , cr in 
km s , plus the fully corrected magnesium index Mg 2 for nearly 
all of them. Distances were taken from the NED database (H = 73 
km/s/Mpc) to compute \ogR e in kpc. We performed analyses both 
with and without Mg 2 (hereafter 4- and 3-parameter configurations 
respectively). We use this parameter as a tracer of the abundance 
of light elements. While it increases as nucleosynthesis proceeds, 
it may not evolve at the same pace as the metallicity (measured by 
Fe/H), because the light elements and Fe are produced in different 
kinds of supernovae arising in stars of different lifetimes. However 
we are mainly interested in an indicator of chemical evolution, not 
in the detailed evolution of elements. No error bars were used in 
this analysis. 



2.2 Cluster analysis 

Cluster analyses do not accept missing values, so that they were 
performed with 696 galaxies in the 3-parameter case (3 galaxies 
did not have redshift hence distance available), and 528 in the 4- 
parameter case (171 galaxies without Mg 2 ). The data for each pa- 
rameter have been standardized to give them uniform weightage as 
we want to classify the galaxies considering all the parameters as 

equally important. 

With the K- means algorithms JHartiganll 19751) . we have ob- 
tained the optimum number of clusters given by the value of K 
associated with the largest jump and partitioned the galaxies ac- 
cording to this grouping (see Appendix lAlfor more details). 



2.3 Cladistic analysis 

We used the 699 galaxies for the cladistic analysis because for miss- 
ing values the algorithm simply guesses values that yield the most 
parsimonious diversification scenario. In this way, cladistic analy- 
sis also provides predictions for undocumented variables. The stan- 
dardization is done automatically through the discretization of the 
variables (Appendix [B}. 

Even though the four parameters show mutual observational 
correlations, they characterize a priori distinct and independent 
physical properties that are expected to evolve and for which some 
sort of evolutionary states could be defined. Note that this holds 
even when the variables are physically linked by the virial the- 
orem, because the virial equilibrium is a physical state that can 
evolve. These parameters are thus adequate for a cladistic analy- 
sis, in which case they are called "characters". 

The cladistic analysis yields trees from which groups can be 
defined. Groups on a cladogram are theoretically defined as evo- 
lutionary groups (or "clades") that comprise a node and all its de- 
scending branches (see Sect. l3.lt . The identification of groups de- 
pends on the rooting of the tree, and must be analysed in regard 
with properties and other characteristics. The trees of the present 
study are rooted with objects or group of objects having low Mg 2 . 

Objects cannot reasonably be described as "more evolved" or 
"less evolved" than others because we do not have an evolutionary 
clock. In addition, the measure of evolution is a very difficult con- 
cept in a multivariate space. This is why "diversification" is a more 
appropriate term. A very basic illustration of this point on the hier- 
archic al growth of dark matter halos can be found in lFraix-Burnetl 
(2009). See also Appendix IB] and Sect. |3.2| for some discussion on 
the measure of diversification. 
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cluster analysis K=4 
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cluster analysis K=4 
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cluster analysis K=7 
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Figure 1. Groups obtained from cluster and cladistic analyses. Left column: 
with three parameters: logR,,, < /i e >, logo". Right column: with the same 
three parameters plus Mg 2 . Top and middle row: cluster analysis results 
for the two possible numbers of groups as given by the K-means analysis 
in each case. Bottom row: cladistic results. The colour coding of groups is 
specific to each box. The colours for the cladistic analysis with 4 parame- 
ters (bottom right box) are the same as in Fig. lBll The grey colour of some 
points on the two upper panel to the right corresponds to missing Mg 2 data. 
Group number in the other plots is used in Fig. ff] Group number "0" cor- 
responds to galaxies that have missing Mg 2 values and were excluded from 
the cluster analyses with 4 parameters. 



3 RESULTS 

3.1 Structures in the fundamental plane 

The results of our analyses are synthesized in Fig.Q] where the dis- 
tribution of the groups in the log cr vs < fi e > plane is presented, for 
the 3-parameter case on the left side, and the 4-parameter case on 
the right side. The K-means analysis finds two peaks, correspond- 
ing to the two upper rows of Fig.Q] at K=4 and K=6 with 3 param- 
eters (without Mg 2 ), and at K=4 and K=7 with 4 parameters (with 
Mg 2 ). As described in Appendix[A] in both cases, the K=4 peak is 
lower than for K=6 or K=7. From a purely statistical point of view, 
the two latter numbers of groups (middle row of the figure) are 
thus favoured. The cladistic result is shown on the bottom row of 
Fig.Q] Five groups can be distinguished on the tree obtained in the 
3-parameter case, 4 of them being very probably true evolutionary 
groups, the other one being an ensemble of successive branches. In 
the 4-parameter case, 7 groups are identified (Sect. [3~2"l l. 
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Table 1. Main correspondence, drawn from Fig[2] between groups from cladistics with 4 parameters and the other groupings. In this table, 
contributions smaller than about 15% are ignored. Note that some groups are roughly subdivided, like C6a, C6b. 
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Figure 2. Histograms showing the distribution of the groups from cladistics 
with 4 parameters for each of the analyses. Same arrangement of plots as in 
Figd The colours of the groups are the same as in Fig. IB II 



It is obvious from Fig.Q]that a U our analyses agree to a small 
number of groups (4-7) and yield very similar groupings. In partic- 
ular, all the groups define separate regions on the plane, not across 
its thickness. It is important to keep in mind that all the groupings 
mentioned in this paper have a statistical meaning, that is their bor- 
ders are not deterministic, and they are defined in a multivariate 
space so that overlapping can be important on 2D or 3D projec- 
tions. 

There is a very good similarity between the 3- and 4-parameter 
configurations for a given statistical method. For the cluster anal- 
ysis, with K=4, there is nearly no difference between 3- and 4- 
parameter results. For K=6 and K=7, the inclusion of Mg2 splits 
two of the groups (red and violet-red) into three (red, green and 
violet-red) along log cr and also slightly along < \i e >. The yel- 
low group for K=6 is reduced in K=7 (cyan) because its low-log <r 
part have many objects with missing Mg 2 indices (represented as 
grey on Fig.[TJ thus excluded from the 4-parameter analysis. For 
the cladistic results, the difference appears to lie along < fi e > in 
the high-log cr part of the plot (one group -purple- split into two - 
magenta and purple- groups) and along log cr for the yellow group 
that is divided into two (yellow and cyan). 



The two top rows of Fig.[T|show that the cluster analysis tends 
to divide the fundamental plane in more or less equal regions, sug- 
gesting a grid which is tilted with respect to the two axes log cr and 
< fi e >. Schematically, the four cases of the cluster analysis show 
the same structure: the high log cr part of the plot is divided in 3 
or 4 groups, mainly but not exactly along < fi e >, while the low 
logo - part is divided between 1 and 3 groups. The dividing bor- 
der between these two parts of the plots is more or less diagonal. 
This behaviour is particularly obvious with K = 4, and the borders 
appear fuzzier in the case with 4-parameter and K=7. 

This trend is still present but less obvious for the 4-parameter 
results (bottom row of Fig. [T}. Noteworthy, adding Mg 2 does not 
reinforce any grouping trend along log cr, that is the apparent grid 
is still tilted in the same direction, which could have occurred in 
the cladistic analysis if the two characters were redundant. Borders 
are somewhat fuzzy, but the most notable difference is that the high 
log cr part of the plot is divided into only 1 or 2 groups. The group 
C5 overlaps several other groups in the middle. These slight differ- 
ences between the two kinds of methods reflect the difference in the 
classification philosophy. But still, the distribution of the groups is 
essentially the same in all cases. In the high log cr part of the plot, 
the 4-parameter cladistic result is thus closer to the cluster analysis 
results than the 3-parameter one. In the low log cr part, the group 
C4 (cyan colour) of the 4-parameter cladistic analysis is nearly the 
same as the cyan group of the 4-parameter cluster analysis with 
K=7. The group C2 (yellow) is composed of galaxies with missing 
Mg2 indices that are not classified in the 4-parameter cluster anal- 
yses. From this point of view, the 4-parameter cladistic and cluster 
results are also quite similar. 

A slightly more precise comparison is made in TableQ]for the 
principal correspondences between groups, and a more quantita- 
tively one is presented in Fig. [2] It must be noted that the group C5 
is often split in two parts that match nearly perfectly the two small 
subgroups visible in Fig. lBlK see Sect. l3.2t . 

Such an agreement between results obtained with two different 
techniques and two sets of parameters gives a high level of confi- 
dence in the structures of the FP of early-type galaxies that we find. 
Since from a statistical point of view, the K=4 cases are less sup- 
ported in the cluster analyses, we will not discuss them any further. 
Since the 4-parameter cladistic result is slightly more robust than 
the 3-parameter cladistic analysis (Appendix [B}, and it is also more 
similar to the cluster analysis results especially with 4 parameters 
and K=7, we concentrate on the 4-parameter cases in the rest of 
this paper. Finally, since the cladistic analysis additionally yields 
the evolutionary relationships between groups, we will mainly dis- 
cuss these 7 groups C1-C7 (Fig. lBlt in the following. Nevertheless, 
to be complete and allow the reader to check that our interpretation 
is not limited to one specific analysis, but relevant for all of them, 
we present diagrams for the other groupings in Appendix Icl 
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Figure 4. Three-dimensional views of the individual fundamental planes. The grey plane on the two bottom right plots is the global fundamental plane. The 
colours of the groups are the same as in Fig. lBll 
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Figure 3. Projection of the tree in Fig. IB II on the log <x-/i e plane. Each 
small branch leads to the group average given in Table|5] The colours of the 
groups are the same as in Fig. lBll 



3.2 Evolution within the fundamental plane 

The groups from cladistics are defined from the most parsimonious 
tree found by the analysis. These groups are only a part of the in- 
formation provided by the tree, that is shown in Fig [BT] for the 4- 
parameter case. It depicts the evolutionary relationships between 
the groups, and has been rooted with CI to orientate the direc- 
tion of diversification, that is the evolutionary distance from more 
primeval classes of objects (see Sect.|2j. On the tree, at least 8 sets 
of branches can be seen, but only 7 are identified as groups for this 
paper (CI to C7). 

As shown by the colours of the branches in Fig. lBll CI has 
the lowest mean Mg 2 . It has been chosen to root the tree (Sect. [2}. 
Following the course of diversification, a split occurs with a branch 
leading to C2 which is slightly more metallic. The node at the split 
indicates a divergence in the evolution of some or all of the proper- 
ties: C2 shares a common ancestorship with C3 to C7 - low metal- 
licity progenitors that could resemble CI objects - but has unique 
characteristics that we must identify in the multivariate space of 
evolutionary parameters. C3 to C7 are more diversified than C2 
with respect to CI. 

Very similar to C 1 , group C2 has a relatively low Mg 2 on av- 
erage. These two groups seem to differ mainly in < fi e >. They both 
lie in the region where the global FP is most distorted and dispersed 
(Sect |3.3b - Their FPs are very different from each other and from 
the other ones, both in orientation and thickness. 

Next in the diversification, the small group C3 is very interest- 
ing because it appears rather early in the diversification scenario but 
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is very similar to groups C6 and C7 in many respects (see Sect. 13.31 
SectEU TableH. 

Group C4 is a well individualized group that is very robust 
since it appears nearly identical in all analyses (Fig.O. It is grouped 
with C2 in two of the 3-parameter analyses, while on the tree, it is 
well separated from the latter, with C3 in-between. The discrimi- 
nation between C2 and C4 thus appears when adding Mg 2 and is 
explained by the evolutionary relation or correlation of this param- 
eter with the three other ones. Note that C2 has a large number of 
undocumented Mgl values (29 out of 70, Tabled, which were not 
taken into account in the 4-parameter cluster analyses. This might 
explain why only C4 appears in all the analyses (Fig.QJ. 

Group C5 is not defined as an evolutionary group but by an 
ensemble of individual branches-galaxies and two small groups. 
Strictly speaking, group C5 is thus not an evolutionary group like 
the other ones. We decided not to individualize these two small 
subgroups to avoid too much complexity in the diagrams, to stick 
as close as possible to the results of the 4-parameter cluster analyses 
that find 7 groups, and because small subgroups could be hard to 
characterize due to a lower statistical significance. Even though we 
have noticed that these two subgroups can be clearly distinguished 
in some of the diagrams, their existence does not modify the main 
results of this paper. It is interesting to note that the two subgroups 
of C5 show up in Fig.[2]and Table[T|as C5a and C5b. The first one 
is next to C3 in the low log cr part of the plot in Fig. Q] and the 
second one is within C6 and C7 in the high log a part. This shows 
that we could have defined 3 groups in this region instead of 2, to 
be compared with the 4 groups of the 4-parameter cluster analysis 
with K=7. 

As already mentioned, C6 and C7 are the two most diversi- 
fied groups in our sample, but their respective places on the tree 
of Fig. lBll could be inverted without modifying any of our conclu- 
sions. They are simply two different groups, more similar to each 
other than to other groups, and that are located at the end of our 
diversification scenario. There is also a subgroup in C7, but we do 
not think it deserves a particular identification in the present paper. 

The evolutionary scenario depicted by the tree in Fig IB II can 
be projected on the global FP (Fig. [3j- The thick track in Fig. [3] 
represents the projection of the tree on the log a -fi e plane, and the 
thin lines represent the beginning of the groups on the tree that are 
here manually extended up to the group average values (Table [3]l. 
Note that distances apparent on this figure are only truncated mea- 
sures of diversification which should be defined in the 4-parameter 
space. In this projection, the diversification within the global FP 
occurs with a main trend along increasing log a and some signifi- 
cant splits along < \i e >. The groups are roughly situated on each 
side of this trend, mainly along < fj e > but also log R e as shown 
on Fig. \5\ C6 and C7 are clearly the two most diversified groups, 
making a radical split, mainly in < fj e >, at the end of the main 
diversification trend. 

The diversification within each groups in Fig.[3]is rather com- 
plex, and we do not describe it in detail in this paper. Globally, 
diversification is always in a direction away from the main thick 
track. 



3.3 Distinct fundamental planes for groups 

In Table|2]we report the coefficients found from a linear regression 
log R e = a log a + b < /i e > + c done within each group and with 
the whole sample. The slope a of the global FP is identical to that 
for group C2, compatible for groups C3, C4 and C6, significantly 
different for C7, and very different for CI and C5. The slope b for 



Table 2. Fundamental planes for the 7 groups and the whole sam- 
ple. obtained by linear regression. The equation of the fundamental plane 
is expressed as log R e = a log cr + b < fi e > + c. The last column gives the 
dispersion about the planes via the standard deviation of the perpendicular 
distance from the corresponding plane. 









a 




b 






std 


CI 




0.67 


±0.13 


0.225 


± 0.032 


-5.4 


±0.8 


0.099 


C2 




1.02 


±0.14 


0.037 


±0.051 


-2.5 


±1.0 


0.070 


C3 




1.32 


± 0.25 


0.349 


±0.071 


-9.3 


±1.2 


0.056 


C4 




1.30 


±0.16 


0.228 


± 0.037 


-6.9 


±0.7 


0.066 


C5 




0.85 


±0.16 


0.305 


±0.019 


-7.3 


±0.6 


0.064 


C6 




1.31 


±0.09 


0.338 


± 0.009 


-9.0 


±0.3 


0.047 


C7 




0.98 


±0.08 


0.349 


±0.014 


-8.5 


±0.2 


0.057 


Full 


sample 


1.13 


± 0.03 


0.334 


± 0.005 
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±0.1 


0.065 
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Figure 5. Statistics of 6 parameters for the groups from cladistics. Each 
boxplot for each group represents the 5 quartiles: median (thick bar at the 
middle of the box), first and third quartiles (limits of the box), minimum 
and maximum (tick bars at the ends of the axis). Outliers are represented by 
circles. The colours of the groups are the same as in Fig. lBll 



the global FP is identical to that for C6, barely compatible for C3, 
C5 and C7, and very different for CI and C4. Group C2 is clearly 
peculiar (it seems to be independent of < /j e >). In summary, only 
C3 and C6 are very close to the global FP, the other ones differing 
more or less, either in one or the other projection or both. C3 and C6 
are only very slightly closer to the virial plane (log R e = 2 log a + 
0.4 < fi e >) but not significantly. 

The dispersion about each plane - the thickness of the plane - 
is similar for the global FP, C2, C4 and C5, but much lower for C6, 
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Figure 6. Histogram of group membership for each galaxy cluster, 
that the name of the cluster appears above each histogram. 



Note 



lower for C3 and C7, and very high for CI (Table [2j. The groups 
with tighter fundamental planes are not necessarily the most diver- 
sified since C3 is one of them. The two least diversified groups, CI 
and C2, are well dispersed and certainly account to a large extent 
for the dispersion about the global FP. 

Fig. [4] presents a 3D view of the planes. Clearly the two less 
diversified groups C 1 and C2, depart the most from the global FP. 
Groups C4 and C5 are also quite different. The two most diversified 
groups, C6 and C7, together with C3, have planes quite similar to 
the global FP. It is however difficult to see a progressive evolution 
in the respective orientations of the planes. We interpret this result 
as an indication that the global FP is not a bent surface, but made 
of a collection of several groups characterizing several fundamental 
planes with different thicknesses and orientations in the parameter 
space. 

3.4 Property distribution of the cladistic groups 

It is important now to understand physically why the groups have 
been individualized by the multivariate algorithms. We again em- 
phasize the fact that these groups have been identified in a mul- 
tivariate space, and that it would be hopeless to try characteriz- 
ing them with one or two obvious properties. It is also important 
to remember that the classification comes from 4 independent, in- 
trinsic and evolutionary pertinent parameters, so that it can be en- 
lightened by various properties even totally different from the one 
used in the multivariate analyses (see our work on galact ic globu- 
lar clusters for a good example, iFraix-Burnet et alj|2009f) . In par- 
ticular, we have computed two quantities that are often used in 
studies of the FP: the dynamical mass M dy „ and the M dy „/L R ra- 
tio by using these two relations: < fj e > = -2.5 log(L R /27rJ?3) + 4.45 
and M dyn a Acr 2 R e /G with A = 3.8 faopkins "et a"Hl2008l) . It is 
important to note that the constant A is empirical and statistical 
for a given sample and is by no way universal. We then obtain: 
\og(M dy „/L R ) =2 2 logcr - \ogR e + 0.4 < fi e > + C, which is the 
"virial" plane if log(M dy „/L R ) and C are constant. The parameter 
C depends on A and some hidden physics, like dark matter, and 
is still not measurable on individual galaxies. We also computed 
M Rabs = -2.5 log(L QR ) + 4.45 =< yu c > -2.5 log(2xrj) + 4.45. Since 
we did not find B and R photometry for the full sample, we did not 
convert these quantities for the B band. 

Average properties are given in Table[3] but boxplots are more 
convenient to visualize the distribution of some property within 
each group and evaluating differences between groups. Fig. [5] 
shows that globally, groups have significantly different properties 
except for redshift and distance to cluster centre, the latter being 
very uncertain due to projection effects. 

A trend toward increase with diversification is present for 
logcr, Mg 2 , and more marginally for M dyn and M dy „/L R . A trend 
toward decrease with diversification seems to be present for < /j e > 
and marginal for M Rabs . There is no trend for log^ c . Note that the 
choice of CI to root the tree of Fig. [BT] implies that the first group 
has the lowest Mg2 on average, but it does not explain the regular 
trends because the relative positions of the groups are independent 
on this choice. 

Galaxies of group CI have the lowest Mgi, they are also large 
(log^ c ) relative to their mass M dy „. They have a low surface bright- 
ness (high < fi e >). CI and C2 show a low logcr, but C2 galaxies 
are smaller, lighter and brighter. 

Galaxies of groups C3 and C6 are on average the largest, the 
most massive and the most luminous. These two groups however 
differ in Mg2, and somewhat in logcr and possibly in the distance 
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to the cluster centre. They have very thin FPs, close to the global 
FP (Sect. [33). 

C4 galaxies have about the same mass and luminosity (M Ra t s ) 
as the ones in CI or C2, but they have a slightly higher logo" and 
surface brightness (lower < /i e >). C5 properties always lie exactly 
between C4 and C6. 

Group C7, which also has a thin FP, is very different from 
C3 and C6. Its galaxies are as small as those of C4, with masses 
equivalent to the average mass of the whole sample. They have a 
high surface brightness (low < fi e >). 

We find that all groups have very similar distribution of T- 
type morphologies (from -5 to -2, not shown), except C6 that has a 
strong excess of T = — 5 galaxies. 

As mentioned above, there is no relation between groups and 
redshift. There is no relation with galaxy cluster either: most clus- 
ters have objects from all groups (Fig. [6^. All have galaxies be- 
longing to either C6 or C7 or both, with a tendency for the largest 
clusters to have more diversified objects. 14 clusters have no galax- 
ies from CI nor C2, but these are clusters with a small numbers 
of objects (less than 8 galaxies). Since we also found very similar 
results with subsamples composed of three well-populated galaxy 
clusters, we conclude that our classification is identical for all clus- 
ters. 

We have performed a Kolmogorov-Smirnov test to check 
whether the distribution of group membership of the galaxies varies 
from cluster to cluster. Among the C? 6 = 1540 possible pairs, only 
66 (4%) have a p- value lower than 0.05 (being equal to 0.012 for 
all the 66 pairs), indicating that there is evidence for a different 
distribution of membership in only 4% of the cases. Interestingly, 
the 66 pairs all have one of the three clusters A0426, A0539 and 
A1656. The different distributions could indicate different environ- 
ments and evolutionary histories for these clusters. 

It appears that the diversification scenario cannot be described 
with only one or even two parameters. Our choice to root the tree 
with Mg2 does not yield a perfectly regular increase of this param- 
eter, and groups that are close in the diversification, like C6 and C7, 
show significant differences in most of the properties. 

3.5 Property correlations of the cladistic groups 

Two-dimensional scatterplots are useful to search for correlations, 
event though they are always projections of a multivariate space 
that naturally induce an apparent dispersion. We show several scat- 
terplots in Fig. [7] and give the correlation coefficients in Table [4] 
for each group. It appears that correlations generally differ between 
groups and the whole sample. 

Some correlations exist for both the whole sample and the 
groups individually. This is the case for the Kormendy relation be- 
tween log^j, and < \i e > that shows an important dispersion for 
the whole sample. However, it is particularly tight for groups C6 
and C7 and rather weak for CI and C4. There is no dependence of 
log R e on < /j e > for C2. It can be concluded that this correlation is 
not universal and only holds for well diversified galaxies. 

The relation between M dy „ and log a is also strongly depen- 
dent on the group (Fig.|7]and Table[4]l. The correlation is generally 
very tight, significantly more than for the whole sample, except 
for C5 and C6 which show an important dispersion. It is impor- 
tant to recall that our computation of M dyn makes it essentially 
proportional to cr 2 R e , the coefficient possibly varying somewhat 
from galaxy to ga laxy or from groups to groups JBolton et al . 2007; 
iHopkins et ai]|2008h . 

Other correlations disappear within each group. The long- 
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Figure 7. Scatterplots discussed in the text. The colours of the groups are 
the same as in Fig. IB II Note that the values of Mgi for the C2 group (yel- 
low) are predictions from the cladistic analysis (see Sect.fS). In the bottom 
right diagram, the full line is a linear fit (slope of 0.51) of th e whole sample 
and the dashed line has a slope of 0.45/1.17 as found by [Robertson et alj 
12006) for pure disk merger remnants (see text for details). 



known global correlation between Mgo and log a appears with a 
correlation coefficient of 0.59 (Table |4j with the sample of the 
present paper. This correlation is similar in C7 but much weaker 
in CI, C2 and C6. It even totally disappears for C3, C4 and C5. 
However, all groups are aligned according to our evolutionary sce- 
nario of Fig. IB 1 1 the group C6 appearing before C7 (as already 
mentioned, C6 and C7 can be inverted in Fig. IB II without modify- 
ing any conclusion). This result shows that the correlation between 
Mg 2 and log <r is historical or spurious and not physical: it is the 
diversification itself that follows this track. 
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Table 3. Average properties of the seven groups of galaxies defined in Fig. IB II Ngal is the number of galaxies per group, 
NoMg2 is the number of missing Mg2 value in the data. 





Ngal 


NoMg2 


log cr 


<fle > 


log« f 


Mg 2 


log(M rf) „) 


M d y,JL R 


M Rabs 


CI 


80 


31 


2.047 ±0.1 14 


20.39 ± 0.45 


0.58 ±0.16 


0.235 ± 0.024 


10.63 ±0.31 


90 ± 38 


-17.49 ±0.71 


C2 


70 


29 


2.054 ± 0.091 


19.56 ±0.24 


0.33 ±0.14 


0.259 ± 0.020 


10.39 ± 0.29 


75 ±27 


-17.06 ±0.69 


C3 


39 


10 


2.284 ± 0.074 


20.10 ± 0.26 


0.78 ±0.19 


0.271 ±0.014 


11.31 ±0.31 


120 ± 32 


-18.78 ±0.77 


C4 


110 


21 


2.173 ±0.072 


19.05 ± 0.30 


0.28 ±0.17 


0.265 ± 0.026 


10.57 ±0.29 


89 ±29 


-17.32 ±0.73 


cs 


111 


24 


2.309 ± 0.054 


19.51 ±0.46 


0.60 ±0.16 


0.288 ±0.013 


11.16 ± 0.19 


119 ±29 


-18.47 ±0.51 


C6 


140 


24 


2.393 ± 0.079 


19.99 ± 0.72 


0.87 ± 0.27 


0.316 ±0.017 


11.60 ±0.35 


147 ± 35 


-19.34 ±0.79 


C7 


149 


32 


2.342 ± 0.091 


18.64 ±0.54 


0.33 ±0.25 


0.301 ± 0.020 


10.96 ± 0.40 


119 ±36 


-18.00 ±0.83 



Table 4. Correlation coefficients for Fig. [7] p- values are indicated in parentheses only when it is higher than 0.001. 





all 


CI 




C2 


C3 


C4 


C5 


C6 


C7 


log R e vs < /i e > 


0.59 


0.24 




0.03 (0.2) 


0.55 


0.34 


0.63 


0.79 


0.80 


logo-vslog(M^„) 


0.72 


0.75 




0.88 


0.85 


0.81 


0.31 


0.47 


0.78 


Mg2 vs log <T 


0.59 


0.23 




0.24 (0.001) 


0.01 (0.6) 


0.008 (0.4) 


0.08 (0.008) 


0.23 


0.53 


M Rabs vs \og{M dyn ) 


0.91 


0.64 




0.69 


0.87 


0.79 


0.73 


0.92 


0.90 


M dyn /L R vsM Rabs 


0.16 


0.08 


(0.01) 


0.09 (0.01) 


0.007 (0.6) 


0.08 (0.003) 


0.15 


0.05 (0.009) 


0.11 


M dyn /L R vslog(M dvn ) 


0.42 


0.11 


(0.003) 


0.07 (0.03) 


0.08 (0.09) 


0.03 (0.07) 


0.02(0.1) 


0.24 


0.37 


logo-vsM Rabs 


0.54 


0.26 




0.39 


0.57 


0.44 


0.08 (0.002) 


0.36 


0.56 


\ogR e vsiog(M dyn ) 


0.73 


0.53 




0.79 


0.90 


0.87 
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0.85 


0.89 
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Figure 8. Diagram showing log R e as a function of 2 * log <x + 0.4 * /i f often 
used to represent the tilt of the FP with respect to the virial plane (dashed 
line). The dotted line is the fit for C2 which has the lowest slope of all 
groups. The slope tends to increase with diversification (see Table|5J. 



This is also true for the correlation between log(M rfv „) and 
M Rabs which is very tight, with similar correlation coefficients for 
the whole sample, C3, C6 and C7. However, the correlation is also 
largely defined by the succeeding groups along the correlation, C6 
being the farthest at one end, and C2 and C4 at the other end, with 
C3 is between C7 and C6. 

Similarly, M dy „/L R is roughly correlated to M Rabs and M dy „, 
but more importantly the relation between these two parameters 



strongly depends on the group. Indeed the first correlation is al- 
ways weak, being slightly higher for the whole sample, C5 and C7. 
The groups seem to be arranged along the trend according to our 
diversification scenario, this trend delineating the global correla- 
tion. The second correlation between M dy „/L R and M dy „ is not very 
strong but still much tighter for the whole sample than for individ- 
ual groups. The succession of the groups is present as well. 



The Faber- Jackson correlation between M Rabs and logo" is a 
mixed case. It clearly depends on the group, being similar for the 
whole sample, C3 and C7, but much weaker for the other groups. 
This correlation also roughly holds for the centroids of the groups 
which appear aligned along the global correlation according to their 
diversification rank. 



A well known global correlation is found between logR e and 
M dyn , with a high significance (Table [4) and a slope of 0.51 (Ta- 
ble[5] see Sect. |4]l. If we now look at individual groups, it is strik- 
ing to see that they all follow a tighter linear relation, except for C 1 
that is slightly more dispersed. The correlation coefficient for C5 is 
equivalent to that of the whole sample and very high for C3, C4, C6 
and C7. Even more remarkably, they are very well separated from 
each other (Fig. [7J- Groups CI, C2, C4 and C7 are stacked per- 
pendicularly to the correlation in an order that closely follows the 
diversification scenario, going basically from the low masses and 
high radii toward higher masses and lower radii (see also Table |5)- 
Groups C3 and C6 are displaced toward the upper right of the cor- 
relation (high mass and high radius), and C5 is in an intermediate 
position. Each group occupies a finite region in the plot, both along 
and across the correlation. Since this result is important, we show 
in Appendix [C] that it is true for all groupings found in the present 
work. Finally, the slopes of the individual correlations (Table[5} are 
lower for C 1 and C2, and then clearly increase with diversification 
being always larger than the global slope. 
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Table 5. log R e as a function of 2 * log <x + 0.4 < /i e > and log Mj y „ . 
The second column gives the slope of the linear regression fit to log R c , vs 
2 * log <x + 0.4 < /i c > which measures the tilt of the fundamental plane with 
respect to the virial plane. The third column gives the slope a of the corre- 
lation: logS f vs logMrf,,,, shown in Fig. [7] Coefficients for dissipationless 
and dissipational systems comes from Robertson et al. 12006) corrected for 



- 0.053). This also implies that our classification, not the FPs, is 
universal within these limits. 



M dy , 



I . The two last columns give the median values for log R e and 



log Mj y „ in the corresponding group. 



Group tilt / virial slope a median median 

log(Re) log(M dyn ) 



CI 


0.42 


0.38 


0.55 


10.65 


C2 


0.40 


0.42 


0.37 


10.47 


C3 


0.74 


0.56 


0.74 


11.27 


C4 


0.61 


0.55 


0.29 


10.61 


C5 


0.71 


0.69 


0.59 


11.17 


C6 


0.81 


0.71 


0.85 


11.60 


C7 


0.70 


0.60 


0.36 


10.98 


Full sample 


0.70 


0.51 


0.51 


10.99 


dissipationless 


0.96 


0.38 






dissipational 


0.83 


^0.49 







4 DISCUSSION 

4.1 Structures in the fundamental plane space 

The first important result of this paper is that we find structures 
in the global FP of this sample despite that it was established to 
be very homogeneous. FP properties are known to differ from one 
sample to another, depending on kinds of objects and redshifts. 
They mostly differ in relative shift and somewhat in inclination. 
This implies that the global FP does not relate galaxies of the same 
kind, but rather traces a diversification scheme of various groups 
of homologous galaxies. A priori, the classification we find is by 
no means universal. It is only valid for the sample used described 
with the 4 parameters \ogR e , < \i e >, logo" and Mgi. Other anal- 
yses have to be performed on other samples and results compared. 
We must mention that we find the same kind of result with a totally 
distinct sample of 500 nearby early-type galaxies (Fraix-Burnet et 

al., i n prep). 

iD'Onofrio et al.l (2008) have performed an analysis of the 
FP for 59 galaxy clusters with redshifts between 0.04 and 0.07, 
quite similarly to our study (4 clusters are in common: A0548SW, 
A2657, A3558, A3716). They compare the FPs obtained for each 
cluster and find that they are not compatible with a universal FP. 
However, they find that the FP properties depend strongly on the 
sample analysed. They use different subsamples, and find that prob- 
ably the different ways early-type galaxies were selected have a 
prominent influence on the fit of the FP. This could have serious 
consequences on the physic s derived from the FP properties, like 
the structure of dark matter JBorriello et alj|2003f) . In our point of 
view, any a priori selection criterion induces an unavoidable bias, 
even if enough care is taken for homogeneity or completeness. A 
multivariate classification is the best objective tool to correlate clus- 
ters with FP properties and investigate whether the FP can be uni- 
versal. In our approach, the groups we have found solely depend 
on the intrinsic properties of the galaxies, independently of cluster 
or type selection. It appears that they do not correlate with clus- 
ter (Fig. [6} and the distributions are identical in 96% of the cases 
(Sect. 13.4b . This implies that both the populations of galaxies and 
the environments in our clusters do not differ much and the evo- 
lution of galaxy properties, that is the diversification of galaxies, 
is similar for nearly all clusters within our redshift range (0.007 



iGargiulo et al.l ( 120091) argue for a bent plane instead of a con- 
tinuous "distribution" of planes depending on some arbitrary pre- 
selection criterion (size, mass, type, ...). ID'Onofrio et al.l (2008) 
find that the plane is actually a bent surface "approximated by dif- 
ferent planes depending on the different regions of the FP space oc- 
cupied by the galaxies of the samples under analysis.". Our result is 
clearly against a universal bent plane or surface, because even a ho- 
mogeneous sample like the one we used, appears to be divided into 
several planes of slightly different inclinations. More significantly, 
our cluster analysis and cladistic results show that the multivariate 
space is divided into different regions, related by the diversification 
process, in each of which a correlation might or mi ght not exist. 
Somewhat in agreement with D'Onof rio et alj ( 120080 . but for dif- 
ferent reasons, we doubt that finding universal and unbiased coef- 
ficients for the FP is a realistic objective. We conclude that what 
is called the global FP is indeed a collection of regions in a mul- 
tivariate space, with projections on the corresponding 3-D space 
more or less planar, depending on the corresponding group. Each 
group can be considered as truly "homologous', as defined in the 
introduction ("similarity due to same class of progenitor") because 
it results from a multivariate analysis, not from arbitrary selection 
criteria. 

Even when the FP is carefully built, a significant dispersion 
remains that suggests another parameter is at play. For instance, it 
has been found that this dispersion co uld be expla ined by the mass- 
distribution depending on mass (Nigoche-Netro et al. 2009). Even 
though our determination of the mass is based on the virial theorem, 
our result seems in agreement with this interpretation. Fig[5jhows 
that the variance within a group is larger for groups of more mas- 
sive galaxies (C3 and C6). Unfortunately, mass, like most of phys- 
ical parameters, is very difficult to measure (see a discussion in 
iHopkins et alj|2008h . However, we think that the dispersion about 
the FP is naturally due to cosmic variance within homologous pop- 
ulations. In our approach, we conclude that mass-distribution de- 
pends on groups, and not necessarily and only on mass. We note 
that for essentially all the parameters considered in this paper, we 
find that intragroup variance is often larger than intergroup differ- 
ences (e.g. FigO. This is certainly intrinsic to the nature of galax- 
ies, and obviously calls for multivariate classifications. In partic- 
ular, this shows that making binary classifications (like low- and 
high-masses objects) may introduce an important bias physically 
unjustified. 

The variation of the ratio M/L with z (i.e. due to passive evolu- 
tion) is often used to explain offse ts with respect to the FP that ar e 
observed at higher redshifts (e.g. Ivan Dokkum & Stan ford 2003). 
This is different from the tilt which can be explained by a variation 
of M/L with L. But assuming M/L to be constant (with a given tilt) 
for a given redshift implies a very strong physical and evolutive 
assumption. This parameter is computed by using the three param- 
eters of the FP space plus the essential assumption of a constant 
ratio between dynamical mass and real mass for all galaxies. In this 
homogeneous sample used in our study, the ratio M dyn /L R is not 
constant across and within the plane, and is different depending on 
the group (Fig. [5] and Fig.|7J. In addition, we do not find any de- 
pendence of this parameter with z but our range in redshift is small. 
Our wo rk thus implies that M/ L is not constant for z = (Sect. 13. 41 
see also llorgensen et al.lll996l) . 
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Figure 9. Diagram showing our groups in the k-space as defined by 
iBender et alj 119921) . The colours of the groups are the same as in Fig. IB II 



of galaxies almost exclusively SOs (while Ide Riicke et alj 120051 
cite t hem as for "bulges of spiral galaxies") whereas our values 
f from iHudson et alj|200lh are for the entire galaxy. There are no 
spirals in the sample we use, but galaxies from C7 should some- 
how share a similar assembly history that these bulges. Groups 
C6 and C7 could be pa ralleled respectively to the two kinds of 
ellipticals suggested by iKormendv et alj ((2009): i) the giant el- 
lipticals that are essentially non-rotating, anisotropic, triaxial, lit- 
tle flattened (El. 5), have cuspy cores boxy-distorted isophotes; ii) 
normal- and low-luminosity ellipticals that rotate rapidly, are rel- 
atively isotropic, oblate-spheroidal, flattened (E3), coreless, have 
disky-distorted isophotes, most bulges of disk galaxies being like 
low-luminosity ellipticals. They also suggest that the first ones are 
formed from dissipationless (dry) mergers and the second ones by 
dissipative (wet) mergers. However, we do not believe that our di- 
chotomy C6-C7 could be explained by a simple dichotomy in dis- 
sipation because such diversified objects have quite complex a his- 
tory that cannot be summarized by one kind of assembly event (see 
below). 

In particu lar, interactions are a quite fr e quent event that trans- 
form galaxies. lAguerri & Gonzalez-Garcial (2009J) have simulated 
the effect of rapid interactions on the evolutions of galaxies. We de- 
rive from their k2-vs-kl plots that our groups, except for C 1 and C2, 
seem to be robust with respect to these kinds of interaction, that is 
these events make a galaxy move within its group. This could mean 
that these groups cannot be formed by rapid interactions only. We 
note that their simulations do not have objects similar to C3 and C6 
groups, but this might be due to the limited number of configura- 
tions they present. Anyhow, C3 and C6 galaxies are the most mas- 
sive ones, implying very probably mergers and accretion to have 
taken place. 



4.2 Groups and assembly histories: the k-space 

IBender et aljd 19921) introduced the k-space representation intended 
to depict the fundamental plane relation with more physical coor- 
dinates. Our groups are represented in Fig. [9] and are outlined in 
the k2-vs-kl projection. Note that our k2 and £3 axes depend on 
fi e which is in the i^-band, and thus appear shifted as com pared 
to other diagrams made in the fi-band. [de Riicke et alj d2005l) have 
compared several samples of gala xies divided in fine m orphologi- 
cal classes, partly using data from IBender et alj dl992D . and some 
semi-analytical models of galaxy formation and evolution more ori- 
ented toward dwarf galaxies. It is obvious on their lc2-vs-kl dia- 
gram that the distribution of groups in our classification appears in 
rather good agreement with their morphological classification. This 
is not surprising since 2 of our 4 parameters used in the cluster and 
cladistic analyses are structural (log a and log R e ), while a third one 
(ji e ) is somehow correlated to stellar concentration. 

This agreement is interesting because our classification does 
not require very sop histicated image analysis procedures (e.g. 
IKormendv et al.1 120090 . In addition, our classification is based on 
physical observed parameters that can be more easily compared 
to models and numerical simulations than detailed features in the 
images. Our multivariate analyses can also be extended to more 
observables to explore the complexity of galaxy diversification. 

Inter estingly, the c ompa rison of our k2-vs-kl diagram with 
that from IBender et alj dl992D shows that our two groups C6 and 
C7 match rather well the distinct regions of respectively bright el- 
liptical galaxies and bulges (of lenti cular or spiral galaxies). Note 
that the data of IBender et alj d 19921) are for the bulge component 



4.3 Groups and assembly histories: tilt of the FP 

The difference between the FP and the virial plane is often stud- 
ied as a tilt i n the log R e vs 2 * logo - + 0. 4 < fi e > representation . 
For instance. iBovlan-Kolchin et al.ld2005l) ; lRobertson et alj d2006h 
have performed simulations of several scenarios of merging histo- 
ries. They find that gas dissipation during mergers of disk galaxies 
can explain the tilt, even if Robertson et al. (2006) do not exclude 
other kinematical or photometrical effects. 

If we compare the tilt of the FP for our different groups (Ta- 
ble [5} with their predictions, we find that our groups agree better 
with dissipational mergers with a lot of gas (coefficients of about 
0.8) than with dissipationless ones (coefficient close to 1). Groups 
C3, C5, C6 and C7 are the closest to the virial plane, but still far 
from it. Groups CI, C2, C4, which are the furthest from the virial 
plane, could either be mergers of systems with much more gas, or 
else have evolved by slow accretions or in isolation. Since disper- 
sion is significant w ithin all groups and coef ficients are much lower 
than those found bv lRobertson et alj (|2P06), especially for less di- 
versified ones, dissipation and gas fraction alone in these merger 
scenarios cannot explain the differences between our groups. As 
mentioned in Sect. 14.21 CI and C2 may not be the result of merg- 
ers. 

In addition, as we have shown in Fig. [4] the planes are tilted in 
the three directions of the FP space, so that the correlation between 
log R e and 2 * log a + 0.4 < fi e > cannot be suffici ent to understand 
all the physics of the complex histories of galaxies. iRobertson et alj 
(2006) also give the variations of the plane coefficients a and b 
(see our Tabled as a function of gas fraction. It appears that they 
are not monotonic and thus more difficult to interpret. Our coeffi- 
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cient a is always much lower than that found by the simulations 
(in the range 0.67 - 1.3 as compared to 1.55 - 2.07 in Table 3 
of iRobertson et al]|2006h and our coefficient b in the lower range. 
However, they use the surface density of stellar mass, whereas we 
use the surface brightness density. The relation between the two 
could well distort the plane in the three directions, and this com- 
parison is probably flawed. I n addition, it seems t hat the slope a 
varies much in the literature JBernardi et al . 2003) depending on 
the wavelength band and the fitting method. Finally, the most dis- 
crepant groups in a, CI, C2 are probably not the result of mergers 
so cannot be compared to such simulations. 



4.4 Groups and assembly histories: radius vs mass 

The relation between R e and M smr is a good diagnostic tool for 
merger histories (e.g. Boylan-Kolchin et al. 2005; Ro bertson et al.l 
2006; Hop kins et alj|2009h . It seems that Mj y JM„ nr oc Mj tar with 
A =: 0.17, but it d epends on each galaxy so that this relation has a 
lot of scatter (e.g. IRobertson et al.ll20060 . Because of this unknown 
scatter, we here use Mj y „ oc M' sl ". in Fig.[7jand Table|5] 

We find a global correlation with a slope of 0.51, in agree- 
ment with the S PSS fit for early-type galaxies (0.56/1.17=0.48 
IShen et alj2003f) but steeper than that found in the simulation result 
for p ure disk merger remnants (0.45/1.17=0.38. IRobertson et al.l 
2006). Dissipational systems (their figure 4 for merging of gas-rich 
disk galaxies with dark matter halos, star formation and supernova 
feedback) give steeper slopes, similar to those of groups CI, C2 
but still not as steep as for all the othe r groups especially C5 and 
C6. More puzzling is that, according to lRobertson et alj ( 120060 re- 
sult, the bulge tends to diminish the slope, whereas we find steeper 
slopes for all the groups except CI and C2. However the steeper 
slopes might be due to mergers of n on-disky objects o r the result of 
repeated merging of small systems dShen et alj2003l) . 

Apart from the slope of the log R e vs logM rfv „ correlation, 
we see from Rob ertson et al . (2006) that the groups that are below 
the global relation have older progenitors, that is progenitors with 
properties characteristic of higher redshifts. This would imply that 
groups CI and C2 (assuming they are the result of a merger) have 
formed recently, hence they had not enough time to diversify much, 
while C7 is clearly the remnants of higher redshifts progenitors. It 
is interesting to note that C6, which from several characteristics 
could be considered as more diversified than C7, seems to be the 
remnant of slightly more recent progenitors. But since it occupies 
the high mass and high radius region of the plot, it is reasonable to 
think that its galaxies are the result of several mergers, the last ones 
being more recent than for C7 objects. Similarly, C3 is above C6, 
showing that its galaxies had a more recent merger. This seems to 
contradict the high resemblance between these two groups with C3 
having a lower Mg2 (Fig.O, but metallicity is also governed by the 
composition of the merging galaxies. 

Mergers with wide or bits are also situated above the global 
logR e vs logM rf ,.„ relation JRobertson et alj|200q) . where CI, C2, 
C3 and partly C6 lie. We could then conclude that groups CI and 
C2 are closer to dissipationless remnants of pure disk wide-orbit 
mergers. Together with our discussion in Sect. 14. 2| and Sect. |4.3| 
this seems to confirm the idea that these objects are the remains 
of interactions or monolithic collapse rather than solely mergers. 
C3 seems to lie above C6, suggesting mergers of wider orbits. This 
apparently contradicts its lower logcr (Fig.O but several different 
kinds of mergers probably occurred, with some particularly violent 
for C6 galaxies. 

Other insights are obtained from lCiotti et alj ( 120071) . Indeed, 



the interpretation of their results should be reconsidered for each 
of our groups instead of the whole sample since the R e vs M plots 
show different slopes depending on the models. For instance, they 
find that 10-equal-mass parabolic mergers yield steeper slopes. An 
important result of their simulations is that parabolic dry mergers 
cannot explain the formation of the scaling laws of early-types, 
but w et mergers do, except maybe for very massive galaxies (see 
also B oylan-Kolchin et alj2 005; Hopkins et al. 2009). Dry mergers 
however do not modify these relations once they exist. This means 
in particular that the FP is formed via wet mergers. Our groups 
C3, C6 and C7 show the thinnest FP and consequently should have 
formed in this manner. Since C3 is less metallic, this merger was 
more ancient. Then subsequent dry mergers have increased its mass 
to the level of C6, whereas C7 galaxies were not subject to impor- 
tant mass increase. We however cannot exclude that C3 might be 
the result of an important monolithic collapse with some subse- 
quent low-metallicity dry mergers. These mergers were probably 
slightly less violent in C3 (having a lower logcr), which could be 
explained by their average larger distance to cluster centres than C6 
(Fig-EJ. 

In real observations, the scatter in the log R e vs log M rfl „ cor- 
relation is quite large (Shenetal. 2003; Franx et al. 2008). In- 
deed, it is is much larger for the whole sample than for each 
evolutionary group (Fig. [7] and Table |3J. The observational scat- 
ter is thus explained e s sentia lly by the "stacked" distribution of 
the groups. iFranx et alj d2008h suggest that the velocity dispersion 
(a 1 oc M dy „/R e ) or the stellar surface density (~ M iy „IR 2 e ) might 
explain the scatter of their correlation with colour. Unfortunately, 
we do not have access to colours for most of the galaxies of our 
sample. We find that \o%R e = a * log M dyn + const, with a given in 
Table [5] depending on the group. This naturally yields the relation 
Mdyn x Re that is different from group to group. Consequently, 
our result shows that the scatter in the log R e vs log M iyn diagram is 
mainly explained by the variation of the relation between mass and 
radius, this relation implying a physical interpretation only valid for 
a particular group. Since \ja ranges from 1.41 to 2.63 and is never 
close to 2 (except for the whole sample!), this seems to indicate that 
the neither velocity dispersion nor stellar surface density entirely 
explain the scatter in the relations between log _/?<,, \o%M iyn and 
colour. Thus our groups are probably not homogeneous in colours. 
We finally note that our Fig. [JJ sho ws th at averaging values of R e 
over bins of M dy „ as in lShen et alj |2003j), mixes different popula- 
tions of galaxies and consequently introduces an artificial scatter 
within each bin. 



4.5 Physical interpretation of the diversification scenario 

Even though the above comparisons are limited by the fact that nu- 
merical simulations could not consider all possible scenarios, they 
show that the groups found by cladistics are evolutionary groups, 
each one gathering objects with the same assembly history. We here 
summarize a possible explanation for the origins of the different 
groups. It must be understood that the assembly and interaction his- 
tories of nearby galaxies are necessarily comple x and comprises of- 
ten sev eral significant transforming events (e.g. lFraix-Burnet et al.l 
l2006d) . 

• C 1 : is chosen as the most primeval group because of its low 
average Mg2. Its galaxies are rather faint, and relatively large for 
their masses, with a low log a. They might be the remains of a sim- 
ple assembly through a monolithic collapse with little dissipation, 
and they were probably perturbed by interactions. 
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• C2 ; its galaxies are less massive and smaller than the ones in 
CI, and they have a slightly higher Mg 2 . They are also somewhat 
brighter. They could be the remains of wind stripping of some kind 
of more diversified objects. 

• C3 : large, massive and bright, they have a low metallicity 
Mg2 with a high velocity dispersion but less than for C6. Although 
this is of low significance, they could be the most distant galaxies of 
cluster centres. They could be the remains of an ancient wet merger 
with subsequent low-metallicity dry mergers. We cannot exclude 
that C3 might be the result of an important monolithic collapse with 
some subsequent low-metallicity dry mergers. 

• C4 : its galaxies look much like the ones in CI and C2, but 
more concentrated, with a metallicity slightly higher than for C2. 
Are they simply galaxies in which stellar formation has been con- 
tinuous, or are they CI galaxies that were initially richer in gas 
that has not been swept like for C2 objects? They could also be the 
remnants of several minor mergers and accretion. 

• C5 : they have intermediate properties between the galaxies of 
C4 and those of C6. This group is indeed composed of two small 
evolutionary groups, one closer to C4 on the FP and the other one 
more or less in-between C6 and C7. Are they intermediate objects? 
We think that they could be objects that cannot be fully described 
by only 4 parameters. Their apparently odd position within the FP 
would the be due to a projection effect of a larger multivariate space 
onto our 4-parameter classification. 

• C6 : its galaxies may be the closest to cluster centres, they 
define the tighter fundamental plane. All the parameters of Fig. \5\ 
are the highest of the sample except for/j f . Most of the galaxies are 
probably large ellipticals usually designated as Es. Their history 
might be complex, with many transformation events, they could 
represent a kind of end state of galaxy diversification. They might 
be the remains of many accretions, minor mergers, together with 
more or less dissipational major mergers. However, they seem to 
be the remains of both wet and dry mergers, the most recent ones 
being of the latter kind. 

• C7 : its galaxies also define a tight FP like C6, but they are the 
smallest of the sample, with the highest surface brightness fi e . They 
are very metallic. They seem to be associated with the remains of 
a dissipative (wet) merger, with very little or no dry mergers. They 
could also have formed through minor mergers and accretions. It is 
striking that they occupy the same region in the kl-k2 space as the 
bulges (of lenticular or spiral galaxies). 



In view of the present results, we conclude that the FP relation 
depends on the group and thus on the histories of galaxies. Obvi- 
ously, the tightness of the relation is different from one group to 
the other being much less for less diversified groups. For more di- 
versified groups, the correlation is tight but still very significantly 
different from the virial plane. In a similar manner as for the bi- 
variate plots of Fig. [7] it is possible that the correlation could be 
historical and not physical. In other words, this planar correlation 
between logR e , logcr and /i e might not be a tilted virial plane due 
to dissipation or a particular behaviour of say M/L, but rather a 
parametric correlation between the evolutions of these parameters. 
If this correlation is tighter in more diversified objects, this is prob- 
ably because they are the result of a mixture of several or many 
transformation events like collapse, interaction and merging, dry 
and wet. To investigate this point further, our analyses should be 
redone with additional measured values of mass, luminosity and 
M/L at least. Then a study of the diversification within each group 
should reveal the combined evolutions of all parameters. 



5 CONCLUSIONS 

In this paper, we have reconsidered the study of the so-called fun- 
damental plane of early-type galaxies. We have used two different 
multivariate clustering tools, cluster analysis and cladistic analysis, 
to explore the 4-parameter space logcr, < fj e >, logR e and Mg 2 . 
With both methods, we used both the three first observables, those 
of the fundamental plane correlation, and the four altogether. 

The sample used in our analysis, taken from H udson et al.l 
u200lh . has 699 objects spread into 56 galaxy clusters. In all the 
analyses, 4 to 7 groups are found, which are consistent from both 
statistical and physical arguments. The very good agreement be- 
tween our different analyses provides a good confidence in the ex- 
istence of structures within the FP. We emphasize that no a priori 
criterion is used to select groups of objects, even in the multivariate 
space. In this paper, we focus on the 7 groups (labelled CI to C7) 
defined by the cladistic analysis because it additionally provides a 
diversification scenario linking the groups. Note that since we are 
in a multivariate space, the wording "diversification" is more ap- 
propriate than "evolution" that applies to a single parameter or that 
has the general meaning of "transformation with time". 

The groups define separate regions on the global fundamental 
plane, not across its thickness. In fact, each group shows its own 
fundamental plane, which is more loosely defined for less diversi- 
fied groups. We conclude that the global FP is not a bent surface, 
but made of a collection of several groups characterizing several 
fundamental planes with different thicknesses and orientations in 
the parameter space. In addition, since all groups are present in all 
galaxy clusters, we conclude that our classification, not the FPs, is 
universal within the redshift range of the sample (0.007 - 0.053). 

By design in the cladistic analysis, each group supposedly 
gathers objects sharing a similar history. They are also related by 
evolutionary relationships that represent a diversification scenario. 
By rooting the scheme with the group of least metallicity, we find 
that the two most diversified groups (C6, C7) have the thinnest FP, 
together with an intermediate one (C3). It probably indicates that 
the level of diversity is linked to the number and the nature of trans- 
formation events like collapse, accretion, interaction and merging, 
dry and wet, and that the fundamental plane is the result of several 
such transforming events. 

Our groups have distinct multivariate properties that can thus 
be interpreted in the light of our current knowledge and under- 
standing of galaxy evolutionary processes. Three groups (CI, C2 
and C4) probably did not form by major mergers and have been 
strongly affected by interaction, some of the gas in C2 objects hav- 
ing possibly been swept out. Three other groups, C3, C6 and C7 
have been formed by dissipative (wet) mergers because they fol- 
low a tight FP relation. In C3, this(ese)) merger(s) must have been 
quite ancient because of the relatively low metallicity of its galax- 
ies. But contrarily to C7, both C3 and C6 have subsequently un- 
dergone dry mergers to increase their masses, more violent in C6 
than in C3. Also, in the k-space, the C7 group clearly occupies the 
region where bulges (of lenticular or spiral galaxies) lie. It could 
have formed through minor mergers and accretions. 

It has been recognized that the properties of the FP depends 
much on the sample. Our approach strongly associates the "sam- 
ple" with the "evolutionary group". Since galaxies are evolutive 
objects, homology, a concept often used for FP studies, can be more 
rigorously defined by 'similarity due to same class of progenitor'. 
Since galaxy histories are so complex, only multivariate studies can 
objectively construct homologous groups. Cladistics is indeed de- 
signed to build homologous groups and provides an evolutionary 
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scenario that relates them. It must be understood that the assembly 
and interaction histories of nearby galaxies are necessarily com- 
lex and com prises often several significant transforming events 
Fraix-Burnet et a l. 2006c, e.g.). The interpretation, based on spe- 
cific assembly histories of galaxies, of our seven groups shows that 
they are truly homologous. Our work also shows that multivariate 
cluster analysis is able to find homologous groups even though it 
cannot predict the evolutionary relationships. 

The properties of the 7 groups clearly reveal that they differ 
in assembly histories. Since they have been obtained directly from 
several observables, the interpretation of the result does not depend 
on any a priori classification. In particular, it partly matches re- 
fined morphological classification because we used 3 parameters 
(among 4) that are linked to structural properties of galaxies. How- 
ever our classification is more easily compared to models and nu- 
merical simulations. Our work can be readily repeated with addi- 
tional observables. 

In addition, the diversification scenario relating these groups 
does not depend on models or numerical simulations. The astro- 
cladistic analysis was based on the assumption that the four param- 
eters logi? e , < fi e >, logo - and Mg 2 are evolutive characters, that 
evolve and can characterize states of evolution of galaxies. The as- 
trophysical interpretation of the diversification scenario and of the 
evolutionary groups demonstrates a posteriori that this assumption 
is correct. This is another proof that cladistics can be applied in as- 
trophysics. Since both cluster and cladistic analyses are only valid 
for the sample in study and the variables used, the present study 
will be extended to other galaxy samples with more parameters. 
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APPENDIX A: CLUSTER ANALYSIS 

Cluster analysis is the art of finding groups in data. Over the last 
forty years different algorithms and computer programs have been 
developed for cluster analysis. The choice of a clustering algorithm 
depends both on the type of data available and on the particular 
purpose. Generally clustering algorithms can be divided into two 
principal types viz. partitioning (K means) and hierarchical meth- 
ods. 

A partitioning method constructs K clusters i.e. it classifies 
the data into K groups which together satisfy the requirement of 
a partition such that each group must contain at least one object 
and each object must belong to exactly one group. So there are at 
most as many groups as there are objects (A" <= n). Two different 
clusters cannot have any object in common and the K groups to- 
gether add up to the full data set. Partitioning methods are applied 
if one wants to classify the objects into K clusters where K is fixed 
( which should be selected optimally). The aim is usually to un- 
cover a str ucture that is alrea dy present in the data. The K- means 
method of lMacOueenl ( 119671) is probably the most widely applied 
partitioning clustering technique. 

Hierarchical algorithms do not construct a single partition 
with K clusters but they deal with all values of K in the same run. 
The partition with K = 1 is a part of the output ( all objects are 
together in the same cluster) and also the situation with K = n 
( each object forms a separate cluster). In between all values of 
K = 1, 2, 3, ...n - 1 are covered in a kind of gradual transition. The 
only difference between K = r and K = r + 1 is that one of the r 
clusters splits in order to obtain r + 1 clusters or two of the (r + 1) 
clusters combined to yield r clusters. Under this method either we 
start with K = n and move hierarchically step by step where at each 
step two clusters are merged depending on similarity until only one 
is left i.e. K = 1 (agglomerative) or the reverse way i.e. start with 
K = 1 and move step by step where at each step one cluster is di- 
vided into two (depending on dissimilarity) until K = n (divisive). 
We feel that for the problem under consideration the partitioning 
method is more applicable because 

(a) A partitioning method tries to select best clustering with K 
groups which is not the goal of hierarchical method. 

(b) A hierarchical method can never repair what was done in pre- 
vious steps. 

(c) Partitioning methods are designed to group items rather than 
variables into a collection of K clusters. 

(d) Since a matrix of distances (similarities) does not have to be 
determined and the basic data do not have to be stored during the 
computer run partitioning methods can be applied to much larger 
data sets. 



method developed by Sugar and James (2003) has been used. This 
procedure is based on rate of distortion theory. By 'distortion' 
we mean a measure of within cluster variation. Let X be a p- 
dimensional random variable (where p components are the parame- 
ters under consideration used for clustering) and C\ , Co, ....Ck be a 
set of K cluster centers. For simplicity, in the present study we have 
considered the mean squared error as the measure of within cluster 

variation which is given by dk=-minc u c 2 , ,c K E[(X-Cx) (X-Cx)] 

where C x is the center closest to X. Using the K means algo- 
rithm (Macqueen 1967) we have first determined the structures 

of the sub populations taking K = 1, 2, 3, etc. For each K 

we estimated the value of distance measure d, = £f , M, where 
Mi = j^jF.^L^x® - xff, xf = I"' galaxy observation for the 
j lh parameter in the;'' 1 cluster (/ = 1, 2, ..., n h j=l,2,....,p) , x = av- 
erage value for the /'' parameter in the;''' cluster, p is the number 
of parameters, n is the number of galaxies and n, is the number of 
galaxies in the i' h cluster (i=l,2,....,K). Here x is the estimate of 
the cluster center closest to the observations in that cluster. Here 
d K can be considered as the estimated minimum achievable distor- 
tion associated with fitting K centers to the data. A natural way of 
choosing the number of clusters is to plot d K versus K and and look 
for the resulting curve (known as the distortion curve). This curve is 
always monotonically decreasing. (Initially we would expect much 
smaller drops for K greater than the true number of clusters because 
past this point ). According to Sugar and James (2003) , the distor- 
tion curve when transformed to an appropriate negative power, will 
exhibit a sharp jump at t he true number of clust ers. In fact they have 
proved (equation (7) of ISugar & Jamell2003h that d K p/2 ~ aK/G 
for K > G and ~ for K < G, where G is the actual number of 
clusters and < a < 1. 

The above relation suggests several possibilities for determining 
the actual number G. In particular one can use the jump method 
which estimates G using argmax K [d K - d K , ], the value of K 
associated with the largest jump in the transformed distribution. It 
also suggests that an appropriate value of y would be p/2. In our 
case as the number of parameters is 3 and 4, one can take y as 1.5 
and 2 respectively. 

We have calculated the jumps in the transformed distortion as 
Jk = (d^ 1 - d^ 2 { ). The optimum number of clusters is the value 
of K associated with the largest jump. The largest jump can be de- 
termined by plotting J K against K and the highest peak will corre- 
spond to the largest jump. In the present situation there are jumps 
at K=4 and K=6 for 3 parameter case and K=4 and K=7 for 4 pa- 
rameter case (Fig. lAlt . The largest jumps are at K=6 and K=7 in 
the former and latter cases respectively. 



APPENDIX B: CLADISTIC ANALYSIS 



For K- means algorithms ( IHartiganlll975h the optimum value 
of K can be obtained in different ways. In the present work the 



Cladistics is the subje ct of numerous books (e.g. lWilev et aljl99ll : 
ISemple & Steell2003f) and the astroc ladistic methodology has b een 
presented in detail in several papers (Fraix-Burnet et al. 2006a.b.c; 
Fraix-Burnet 2009; Fraix-Burnet et al. 2009). We refer the reader 
to these references for a complete description of the method used 
here. For the present analysis, we have built a matrix with val- 
ues of the four parameters for all galaxies. The values for each 
parameter were discretized into 30 equal-width bins represent- 
ing supposedly evolutionary st ates (see e.g. Golobo et al. 2006; 
iThuillard & Fraix-Burnetl I2009T) . We adopted the popular parsi- 
mony criterion, which selects the most parsimonious tree among 
all possible arrangements because it represents the simplest evolu- 
tive scenario for the sample. (Maximum parsimony corresponds to 
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Figure Al. Jump Jx as a function of K for the 3- (top) and 4- (bottom) 
parameter cases. 



the minimum number of changes of states for all parameters that 
occur along the paths between all objects. This number is unique 
to each tree.) Fig IB II is a majority-rule consensus tree reflecting 
the most common features in all equally most parsimonious trees. 
For this paper, the maximum parsimony searches were performed 
using the heuris tic algorithm implemented in the PAUP*4.0blO 
dSwofford 120031) package, with the Multi-Batch Paup Ratchet 
method (http://mathbio.sas.upenn.edu/mbpr ). Heuristic methods do 
not explore the parameter space of all possible tree arrangements, 
which would take a prohibitive computer time for hundreds of ob- 
jects, but try to find the minima. They cannot guarantee finding 
the absolute most parsimonious trees but generally require far less 
computer time while being quite effective. The Ratchet approach 
still improves this efficiency. The results were interpreted with the 
help of the Mesquite software (Maddison & Maddison 2004) and 
the R-package (used for graphics and statistical analyses). In con- 
trast to multivariate distance methods, undocumented values are not 
a problem in cladistic analyses. This is why the galaxies missing 
Mg 2 determination have not been excluded. 

In the two cases of the present study, with and without Mg 2 , 
we find rather robust trees. To assess this robustness, we proceeded 
in the following way. We first selected 3 clusters, A 1656 (Coma 
cluster, 56 galaxies, the largest cluster in this sample), A2199 (38 
galaxies) and A3526 (37 galaxies). Cladistics analyses were run on 




Figure Bl. Cladogram obtained with 4 parameters. It is a majority rule 
consensus of 986 equally parsimonious trees. The branches are coloured 
according to the value of Mg2 (from dark blue to yellow through green as 
Mg2 increases). The outer annulus is formed by the names of the galaxies, 
with colours according to the group (numbered from CI to C7) as detailed 
in TablelCil 

these 3 subsamples separately, and then with the 3 together. The 
four tree structures were compatible, indicating that an underlying 
arrangement did exist. This was confirmed by the analysis done on 
the whole sample that converged quite easily despite the low num- 
ber of characters as compared to the number of objects. It must be 
noted that the arrangement of the objects on a tree is constrained 
by the level of information provided by the parameters. When their 
number is relatively low, the constraints are weaker and the number 
of possible most parsimonious trees is larger, making the conver- 
gence toward compatible trees less probable. In the present case, 
we find that the 3-parameter result is slightly less robust. Neverthe- 
less, all trees with and without Mg 2 are compatible. 

The choice of the root of the tree orientates the evolutionary 
processes and the "rank" of diversification of the groups as seen on 
Fig lBll It indicates the ancestral states of the characters. The trees 
of the present study are rooted with objects or group of objects 
having low Mg 2 because we find it is the only objective criterion 
of ancestrality. However, this choice must imply a consistent evo- 
lutionary scenario for all the other parameters, which is the case 
here. 



APPENDIX C: COMPLEMENTARY INFORMATION 

We present here three figures for the other groupings found by the 
other analyses. Fig. lCll is the equivalent of Fig.|5]but for the cluster 
analysis with 4 parameters and K=7 with colours as in Fig.Q](mid- 
dle row, right column). If ones recall that C6 and C7 are split into 
4 groups (4 to 7 on the figure), and that C2 and C4 are mixed to- 
gether in this cluster analysis, the correspondence between the two 
boxplot figures is striking. 

Fig. lC2l Fig. lC3l and Fig. lC4l are the equivalent of Fig.|7]for all 
the analyses. Fig. lC5l shows the tilt with respect to the virial plane 
like in Fig.[8]for all the groupings. 

Finally, we give in Table IC li the list of all the galaxies used in 
this study distributed in the groups CI to C7. 



Structures in the fundamental plane 17 



rS?*j 



a 



^at 



2 3 4 5 6 7 



let *^ 



2 3 4 5 8 7 



2 3 4 5 6 7 



Mg2 



log(Mdyn) 



T '• 




T 


^^ 


?- 


"i 


T Bg:— T 


s_ 


T B.^" T 


■ 


R 


■? " 


-L 







2 3 4 5 6 7 



2 3 4 5 6 7 




2 3 4 5 6 7 



distance to cluster center (kpc) 



X X * 

■t9 ■ 

:B: Tit 

1* 


i- 



fl 




2 3 4 5 6 7 



2 3 4 5 6 7 



ffiBaBlB 



2 3 4 5 6 7 



Figure CI. Same as Fig. [5]but for groups from cluster analysis with 4 pa- 
rameters and K=7. The colours are the same as in Fig.[T](middle row, right 
column). 
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author. 



18 D. Fraix-Burnet et al 



ML 
Be 




8- 

s ; as 



• :.*i«HSW ■ 




'^*i'.. 






Figure C2. Diagrams as in Fig. [7] but for the cluster analyses with K=4.The colours are the same as in Fig.[T](top row). 



M" 



2$r ■ 



ati£ 



:.?** 



5 . • uigg&i ■ ' 1 





&&<•■ ■ 




£&*• 




Bsv 




-■-'-4'^'' 





Figure C3. Diagrams as in Fig. [7] but for the cluster analyses with K=6 and K=7.The colours are the same as in Fig.[T](middle row). 
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Figure C4. Diagrams as in Fig. [7] for the two cladistic analyses with 3 and 4 parameters. The colours are the same as in Fig.[T](bottom row). 
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Figure C5. Tilt with respect to the virial plane like in Fig.[8]for all the groupings. 
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Table CI. List of galaxies for each group. 



CI 



A2199:B-008 

A0576:SMC-A 

N4661 

A1656:D-204 

7S21:PP-S06 

A2877:D-040 

A0189:SMC-I 

A0262:PP-A05096 

10293 

A0539:D-043 

A0957:D-033 

A3733:SMC-G 



A1060:ST-026 

A1177:SMC-B 

A1656:D-107 

A3558:FCP-31 

S0761:FCP-14 

A2063:D-035 

A2199:B-087 

A2634:D-087 

A4038:D-052 

A0262:B-042 

A0576:SMC-J 

A4038:D-033 



E323-008 

A3733:SMC-K 

N0085A 

E243-052 

N0548 

S0301:D-024 

A0426:PP-P07 

A3381:D-037 

A1016:SMC-F 

A1060:JFK-R245 

E322-100 

S0761:FCP-07 



A3558:FCP-35 

A3574:W-050 

A2052:MKV-13 

A2063:D-034 

A2199:B-034 

A4038:D-083 

A4038:D-039 

U02717 

A1016:SMC-E 

14041 

A4038:D-045 



11618 

PISC:PP-Z10020 

A0194:D-071 

S0301:D-034 

A0539:D-054 

A0957:D-049 

A1016:SMC-G 

A1060:JFK-RMH50 

A3526:D-026 

A3558:FCP-56 

E445-054 



A2052:EFR-B 

A2063:D-033 

14767 

A4038:D-067 

A4038:D-037 

A0576:SMC-B 

A1016:SMC-A 

A1656:D-081 

E445-040 

A2877:D-033 

A0189:SMC-J 



U01269 

S0301:D-022 

A0539:D-064 

A0957:D-035 

A1060:JFK-R225 

A1139:D-016 

A3526:D-047 

A3558:FCP-39 

S0761:FCP-24 

A2063:D-065 

A2199:RS-028 



C2 



7S21:PP-S07 

N0386 

A2877:FCP-24 

H0122:PP-H01051 

A0262:B-038 

S0301:D-031 

N2330 

A1060:JFK-RMH26 

E437-045 

E322-099 



A1656:D-087 

A3558:FCP-15 

E384-029 

S0805:FCP-09 

A4038:D-076 

A0957:D-050 

A2806:SMC-F 

PISC:PP-Z01032 

A2877:D-025 

A0194:D-045 



A0569:SMC-L 

A1060:JFK-RMH35 

A1177:SMC-C 

E322-102 

N4729 

Z160-027 

A1656:D-193 

A3570:SMC-64 

S0753:W-051 

El 04-002 



A2806:SMC-D 

N0398 

11638 

A0194:D-052 

S0301:D-026 

A3389:D-043 

A0957:D-037 

E501-049 

A1257:SMC-GC 

A3526:D-050 



14011 

A3571:SMC-40 

S0753:W-095 

A3716:D-061 

A4038:D-053 

PISC:PP-Z01034 

A2877:D-011 

11648 

U01040 

S0301:D-020 



A1060JFK-R338 

A3526:D-040 

A3526:D-027 

A3526:D-015 

A1656:D-156 

A1656:D-027 

S0753:W-017 

A2199:B-044 

A4038:D-060 

A4038:D-049 



A3389:D-053 

A0999:SMC-C 

A3526:D-059 

A1656:D-135 

A4038:D-066 

S0761:FCP-26 

S0301:D-027 

A0400:D-057 

A3526:D-035 

N4743 



C3 



N0420 

A0426:PP-P08 

A3381:D-112 

A1257:SMC-G 

A3570:SMC-50 

A3558:FCP-17 



A2199:B-015 

A0347:PP-B16 

A3733:SMC-I 

10171 

A0539:D-039 

11116 



A0569:SMC-R 

N4616 

A3571:SMC-10 

A2634:D-031 

A0426:PP-P1 1 

A0262:B-018 



N2235 

J8:PP-J03049 

A0539:D-059 

N2329 

N4929 

A3571:SMC-44 



N5304 

A4038:D-043 

A0539:D-042 

A3716:D-098 

11860 



A0548:D-007 

A0957:D-043 

A3558:FCP-57 

MKW12:FCP-09 

A4038:D-065 



13955 

A0539:D-031 

10310 

A3381:D-025 

A0957:D-044 



C4 



11548 

A2877:D-035 

A0539:D-041 

A1016:SMC-C 

A3526:D-049 

A1656:D-238 

N4876 

A2199:B-048 

A2199:B-073 

U01030 

A0548:D-020 

E437-021 

A3537:SMC-156 

E5 10-054 

A2634:D-093 

N0397 



J8:EFR-H 

A0426:PP-P33 

A1060:JFK-RMH28 

A1314:SMC-D 

A3526:D-036 

A1656:D-153 

14133 

S0753:W-047 

A2199:B-033 

S0805:D-029 

A2634:D-104 

PISC:PP-Z01047 

U00996 

A0539:D-075 

E436-044 

A3526:D-009 



S0753:W-010 

A2199:B-054 

A3656:SMC-X 

N0909 

E488-009 

A1257:SMC-E 

A1656:D-096 

A3716:D-081 

A4038:D-059 

PISC:PP-Z01073 

A0194:D-012 

J8:PP-J01080 

A0539:D-052 

A1139:D-041 

A1367:B-020 

E322-101 



A1736:D-144 

S0761:FCP-05 

A2199:B-047 

S0805:D-021 

A2634:L-B03C 

N0388 

U01003 

A3381:D-067 

A1060:ST-034 

A3526:D-033 

13947 

S0753:W-012 

A2199:B-084 

A3656:SMC-S 

A0539:D-063 

10464 



A3571:SMC-29 

E286-029 

A2806:SMC-E 

A2877:D-021 

11696 

A0400:D-017 

A3389:D-048 

A1177:SMC-H 

E322-075 

E323-009 

A1656:D-207 

A1736:D-039 

A2063:D-046 

A2199:B-021 

E235-039 

A2877:D-048 



A3381:D-033 

A1314:SMC-E 

A1656:D-210 

13957 

A3581:SMC-76 

A2199:B-095 

A4038:D-068 

A0539:D-051 

A0569:SMC-B 

12955 

E384-036 

A2634:B-013 

A2877:D-037 

11680 

A0262:B-019 



A1257:SMC-C 

E322-089 

N4824 

N4906 

A3558:FCP-29 

A2199:B-028 

14748 

A2634:D-079 

N0501 

N4850 

N4882 

A1314:SMC-B 

A0426:PP-P21 

A0957:D-046 

A0426:PP-P20 



Table CI - continued 
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C5 


N0394 


U02673 


A3581:SMC-75 


A0539:D-044 


A2063:D-050 


15354 


A3558:FCP-33 


N0385 


12744 


A2063:D-089 


A3389:D-060 


A3716:D-117 


N0759 


E445-059 


N0712 


N4946 


A3656:SMC-I 


A1016:SMC-B 


N0912 


11858 


A2063:D-073 


A0426:PP-P15 


N5424 


A3744:SMC-I 


N4767 


A2063:D-071 


A0576:SMC-I 


A2199:B-026 


A3381:D-075 


A2199:B-005 


A2634:D-068 


N4908 


A2877:D-045 


N4881 


A3733:SMC-H 


A0957:D-054 


A2634:B-021 


11569 


E325-013 


A0400:D-070 


A3558:FCP-21 


A2634:D-038 


N3841 


A0076:D-016 


J8:EFR-I 


E51 1-021 


A0539:D-062 


A2199:EFR-0 


N0079 


A1656:D-161 


A3526:D-041 


A0539:D-057 


N6146 


A0569:SMC-N 


A3733:SMC-B 


U01837 


A3558:FCP-06 


A2877:D-028 


N4816 


A3656:SMC-P 


E501-003 


A0999:SMC-G 


N1224 


A3581:SMC-78 


J8:PP-J07038 


A1736:D-137 


A2634:D-043 


A1656:D-140 


N0382 


N3308 


A2063:D-059 


N1273 


N5438 


A2657:D-070 


A1656:D-230 


N0560 


A1656:D-206 


A2199:EFR-H 


A3381:D-100 


A3716:D-090 


N0541 


A3571:SMC-164 


A0400:D-089 


A3574:W-024 


A3744:SMC-T 


A0999:SMC-E 


A0189:SMC-C 


11806 


A2052:MKV-60 


A0539:D-016 


N6158 


A2634:D-107 


E322-081 


A2063:D-074 


A3381:D-034 


A2199:B-019 


10458 


A3744:SMC-Q 


A0076:D-018 


N4864 


A2877:D-042 


N4854 


A3716:D-099 


A1139:D-030 


N4683 


J8:EFR-K 


A3558:FCP-04 


11807 


A3558:FCP-07 


A2634:D-130 


A1656:D-239 




A2806:SMC-C 


N6998 


N3309 


L.0 

A2063:D-072 


N2340 


E325-004 


N0410 


A0400:D-058 


15341 


A1314:SMC-A 


A2634:D-074 


N4696 


E5 10-066 


N0564 


A0539:D-045 


A2657:D-031 


A1656:D-240 


A4049:D-055 


E444-046 


N6173 


N1278 


A0957:D-059 


N0212 


E509-008 


11566 


N6999 


A3716:D-080 


A0999:SMC-D 


S0753:W-073 


A0539:D-047 


A3558:FCP-13 


E243-045 


N0703 


E464-018 


12738 


A3716:D-067 


10661 


14329 


N0508 


10313 


N7728 


N3862 


A2634:B-016 


N4889 


E5 11-026 


N1272 


A0569:SMC-Q 


15353 


N4926 


N0083 


N6160 


El 87-020 


N2230 


N4869 


J8:EFR-D 


A3558:FCP-03 


N0379 


E243-049 


A3733:SMC-C 


A1139:D-039 


A2199:B-020 


U05515 


E325-016 


A0189:SMC-A 


10312 


A2634:D-057 


A1314:SMC-G 


15342 


N4839 


14374 


U01308 


A0539:D-050 


A2657:D-071 


14051 


A4049:SMC-E 


U09799 


14765 


A0539:D-068 


12597 


N0708 


A3558:FCP-02 


11568 


15358 


A3716:D-078 


10613 


A2063:D-090 


A3381:D-055 


A3558:FCP-16 


11633 


A0347:PP-B07 


N7016 


10708 


E286-049 


N3551 


N5419 


N0545 


N1293 


N7735 


N4709 


A2657:D-043 


A3558:FCP-05 


E5 11-032 


11907 


A0576:SMC-D 


15362 


N4927 


11565 


N6166 


A3716:D-065 


A0576:SMC-C 


N4923 


J8:EFR-C 


A3558:FCP-18 


N0383 


N0499 


A3733:SMC-A 


10660 


E104-007 


N3311 


A3571:SMC-171 


N0507 


N1283 


A2634:D-077 


N3842 


A2634:D-036 


N4874 


E5 11-023 


11803 


A3381:D-021 


A2657:D-064 


E443-024 


N0080 


A2063:D-060 


14931 


A3381:D-056 


13959 


U01841 

CI 
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S0761:FCP-04 
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A3526:D-046 
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A3558:FCP-14 


A2634:D-102 


A1139:D-029 
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A0539:D-048 
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A1060:JFK-RMH79 


N4919 


A2199:B-061 


E103-046 


N0687 


A3389:D-049 


A1228:SMC-M 


N4730 


A3558:FCP-25 


A3716:D-141 


A3716:D-084 


A0400:D-044 


A1060:JFK-RMH29 


N3851 


N4872 


A3558:FCP-09 


A2634:D-075 


A2634:D-119 


U02698 


E323-005 


N4860 


14021 


A3571:SMC-32 


A2634:D-069 


A4049:SMC-D 


A0548:D-051 


A1656:D-136 


A3558:FCP-34 


S0753:W-049 


N5397 


A4038:D-055 
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A0548:D-017 
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N7014 
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E436-045 
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A2199:B-030 
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E382-002 


A2634:D-080 


A3716:D-051 
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A2657:D-072 
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E243-041 


S0301:D-017 


A1139:D-037 


A3574:W-074 


A3571:SMC-13 


A4038:D-051 
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N0547 


A0426:PP-P26 
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A3581:SMC-77 
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N0392 


A0539:D-049 


N0911 
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10709 
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A2199:B-050 
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E437-013 
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A2634:B-030 


E235-049 
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